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Q_i' Abstract 

We investigate the structure of the D s i (2460, 2536) (J p = 1 + ) mesons via analyz- 
ing the semileptonic B c — > D s \l + l~, I = T,fi,e and B c — > D s \vv transitions in the 
framework of the three-point QCD sum rules. We consider the D s % meson in two 
ways, the pure \cs) state and then as a mixture of two \ 3 Pi) and I 1 -Pi) states. Such 
type rare transitions take place at loop level by electroweak penguin and weak box 
diagrams in the standard model via the flavor changing neutral current transition of 
b — > s. The relevant form factors are calculated taking into account the gluon con- 
densate contributions. These form factors are numerically obtained for \cs) case and 
plotted in terms of the unknown mixing angle 9 S , when the D s \ meson are considered 
as mixture of two | 3 -Pi) and I 1 Pi) states. The obtained results for the form factors 
are used to evaluate the decay rates and branching ratios. Any future experimental 
measurement on these form factors as well as decay rates and branching fractions and 
their comparison with the obtained results in the present work can give considerable 
information about the structure of this meson and the mixing angle 8 S . 
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1 Introduction 



The structure of the even-parity charmed D s j mesons has not known exactly yet and 
has been debated in the quark model. The observation of two narrow resonances with 
charm and strangeness, -D s o(2317) in the invariant mass distribution of D s 7r° [1-6] and 
-D s i(2460) in the D*7r° and Dg'j mass distributions [3-8], has raised discussions about the 
structure of these states and their quark contents [9, 10]. Analysis of the _D So (2317) — > -D*7, 
L> s i(2460) -> £> s *7 and £> Sl (2460) -> D S0 (2317)7 shows that the quark content of these 
mesons are probably cs [11]. Among these mesons, the axial vector charm-strange meson 
D s i is more attractive ones, because the discovery of the _D sl (2460)( J p = 1 + ) meson [2- 
5] and its measured mass indicated a lower mass than expected in potential model (PM) 
[12] and quark model (QM) [13,14] predictions. In other words, the .D s i(2460) does not 
fit easily in to the cs spectroscopy [15]. However, some physicists presumed that this 
discovered state is conventional cs meson [16-26]. Many different theoretical efforts have 
been dedicated to the understanding of this unexpected and surprising disparity between 
theory and experiment [27-35]. As a result of the above discussion, we will consider the 
D s i meson in two ways, the pure \cs) state and also as a mixture of two \ 3 P\) and 1 1 -Pi) 
states. 

Heavy-light mesons are not charge conjugation eigen states and so mixing can occur 
among states with the same J p and different mass that are forbidden for neutral states. 
These occur between states with J = L and S = 1 or [15]. Hence, the mixing of the 
physical D s i and D' sl states can be parameterized in terms of a mixing angle 9 S , as follow: 

\D 8l ) = smds^P^ +cos9 s \ 1 P 1 ), 

\D' al ) = cos6 s \ 3 P l ) -stn6 s \ 1 P l ). (1) 

where, the spectroscopic notation 2S+1 Lj has been used to introduce the mixing states. 
Considering | 3 Pi) = |Z? sl l) and I 1 Pi) = |Z? sl 2) with different masses and decay constants 
[36], we can apply these relations for axial vectors -D s i(2460) and _D sl (2536) mesons with 
two different masses, i.e., 

|D 8 i(2460) ) = sin6 s \D sl l) + cos9 s \ D sl 2) , 

|D s i(2536) ) = cosO s \D sl l) - sin9 s \D sl 2) . (2) 

The masses of D sl l and D s i2 states are presented in Table 1. These values have been 
obtained in QM approach. 



Ref 


[12] 


[13] 


[14] 


£ s il( 3 Pi) 


2.57 


2.55 


2.535 


Ds^P,) 


2.53 


2.55 


2.605 



Table 1: Masses of l 1 P 1 and l 3 Pi heavy-light mesons in quark models. 

Note that, in the heavy quark limit the physical eigen states D s \ and D' sl can be iden- 
tified with P\^ 2 and P^ 2 with notation Lj, where j is the total angular momentum of the 
light quark [37], corresponding to 9 S = —54.7° [36]. 
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In this work, taking into account the gluon condensate corrections, we analyze the rare 
semileptonic B c — > D sl I = T,/j,,e and B c — > D sl vv transitions in three-point QCD 

sum rules (3PSR) approach. Note that, the B c — > (£)*, .D*, D s i(2460))z/z/ transitions have 
been studied in Ref. [38], but assuming the D s \ only as cs. The B c — > D g Z + /~/W [39], 
£? c — > £>* Z + /~, (g = d, s) [40] transitions have also been analyzed in the same framework. 

The heavy B c meson contains two heavy quarks b and c with different charges. This 
meson is similar to the charmonium and bottomonium in the spectroscopy, but in contrast 
to the charmonium and bottomonium, the B c decays only via weak interaction and has a 
long lifetime. The study of the B c transitions are useful for more precise determination of 
the Cabibbo, Kabayashi, Maskawa (CKM) matrix elements in the weak decays. 

The rare semileptonic B c — > D s \l + l~ jvv decays occur at loop level by electroweak 
penguin and weak box diagrams in the standard model (SM) via the flavor changing neutral 
current (FCNC) transition of b — > sl + l~ . The FCNC decays of B c meson are sensitive to 
new physics (NP) contributions to penguin operators. Therefore, the study of such FCNC 
transitions can improve the information about: 

• The CP violation, T violation and polarization asymmetries in b — > s penguin chan- 
nels, that occur in weak interactions, 

• New operators or operators that are subdominant in the SM, 

• Establishing NP and flavor physics beyond the SM. 

To obtain the form factors of the semileptonic B c — > D sl (2460 [2536]) transitions, first, 
we will suppose the _D s i(2460) and D sl (2536) axial vector mesons as the pure \cs) state and 
calculate the related form factors. Second, we will consider the D si meson as a mixture of 
two components \D sl l) and \D s i2) states and calculate the form factors of the B c — > D sl l 
and B c — > D s ±2 transitions. With the help of Eq. (2) and the definition of the form factors 
which will be presented in the next section, we will derive the transition form factors of 
B c — > D sl (2460 [2536]) decays as a function of the mixing angle S . The future experimental 
study of such rare decays and comparison of the results with the predictions of theoretical 
calculations can improve the information about the structure of D sl meson and the mixing 
angle 9 S . 

This paper is organized as follow. In section 2, we calculate the form factors for B c — > 
D s i transition in the 3PSR. In section 3, the two-gluon condensate contributions as non- 
perturbative corrections are calculated. The calculation of the decay rates for B c — > D s \l + l~ 
and B c — > D sl u9 transitions are presented in section 4. Finally, section 5 is devoted to the 
numeric results and discussions. 

2 The form factors of B c — > D s \ transition in 3PSR 

In the standard model, the effective Hamiltonian responsible for the rare semileptonic B c — > 
D sl l + l~ and B c — > D sX vv decays, which are described via b — > s l + l~ loop transitions (see 
Fig. 1) at quark- level, can be written as: 

C 9 e// s 7m (1 - l 5 )bh/ + Cio 57^(1 - i 5 )b£ lfll5 £ 



n 



eff 



2Wi tb ts 



2 



(3) 



where C 7 e// , C 9 e// and C 10 are the Wilson coefficients, is the Fermi constant, a is the 
fine structure constant at the Z mass scale and V^- are the elements of the CKM matrix. 
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Figure 1: The loop diagrams of the semileptonic decay of B c to D sl . The electroweak 
penguin and box diagrams are shown in parts (a) and (b), respectively. 

These loop transitions occur via the intermediate u, c, t quarks. In the SM, the mea- 
surement of the forward-backward asymmetry and invariant dilepton mass distribution in 
b — > q'l + l~ , (q' = d,s) transitions provide information on the short distance contributions 
dominated by the top quark loops [41]. The electroweak penguin involving the contribu- 
tions of photon and Z bosons is shown in Fig. 1(a) and Fig. 1(b) presents the contribution 
of the W box diagram. It is reminded that the b — > s vv transition receives contributions 
only from Z-penguin and box diagrams. 

The transition amplitude of B c — > D sl l + l~/vv decays is obtained sandwiching Eq. (3) 
between the initial and final states, i.e., 



M = -^=V tb V* 



C e 9 ff (D sl (p') I s 7m (1 - l,)b I B c (p))£^£ 



Gpct 

+ C10 (D sl (p') I s 7m (1 - 7s)& I B c {p))l llxl d 

-2 C e 7 ff ^(D sl (p f ) I s ia^il + 75)6 I B c (p))l % 



(4) 



where, p and p' are the momentum of initial and final meson states, respectively, and 
e is the polarization vector of the D s \ meson. Our aim is to parameterized the matrix 
elements appearing in Eq. (4) in terms of the transition form factors considering the 
Lorentz invariance and parity considerations. 

9 aB c ^D s1 , 2\ 

(D sl (p',e)\s W b\B c (p)) = v - {q } e, va ^p a p ffS , 

m Bc + m Dal 

(D sl (p', e) I s llx b I B c (p)) = - iA B ^(q 2 )(m Bc + m Dsl )< + i ^^^ {e*p)P» 

m B c + m D 3l 
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A B^D S1{2) 

+ i — — — £ P <fc , 

"Be + ™D sl 

(£> a i(?/,e) | s^g-'Tsft | S c (p)) = 2 T Bc ~* Dsl (q 2 ) ie^e*W, 

(D sl (p f , e) | 3(7^6 | B e (p)) = T B ^(q 2 )[s;(m 2 Bc - m 2 D J - (s*p)P fl ] 

+ T B ^( q 2 )(e*p)[q,- q \ 2 P M ], (5) 

where Af c ~* Dsl (g 2 ), i = V, 0, 1, 2 and T B °^ sl (q 2 ), j — V, 0, 1 are the transition form factors, 
^/x = (p + P% an d % = (.V ~ p')n- Here, g 2 is the momentum transfer squared of the Z 
boson (photon). In order to our calculations be simple, the following redefinitions of the 
transition form factors are considered : 

A^ D ^(q 2 )(m Bc + m Dal ), 

A B ^(q 2 ) 
m Bc + m Dal ' 

- T B ^(q 2 )(m 2 Bc -m 2 Dsl ) : 

(6) 

To calculate the form factors within three-point QCD sum rules method, the following 
three-point correlation functions are used: 

n^pVV) = ? j ^V^e^(0 I T[J^(y)jV- A (0)J B <\x)] | 0), 

n^V,^ 2 ) = ^/^ye-^^oiriJ^^jJ-^J^^llo), (7) 

where J® sl (y) = 07^753 and J Bc (x) = C756 are the interpolating currents of the initial and 
final meson states, respectively. J^~ A = 57^(1 — 75)6 and Jj~ PT = sa^g^l + 75)6 
are the vector-axial vector and tensor-pseudo tensor parts of the transition currents. In 
QCD sum rules approach, we can obtain the correlation function of Eq. (7) in two sides. 
The phenomenological or physical part is calculated saturating the correlator by a tower of 
hadrons with the same quantum numbers as interpolating currents. The QCD or theoretical 
part, on the other side, is obtained in terms of the quarks and gluons interacting in the 
QCD vacuum. To drive the phenomenological part of the correlators given in Eq. (7), 
two complete sets of intermediate states with the same quantum numbers as the currents 
J Dsl and J Bc are inserted. This procedure leads to the following representations of the 
above-mentioned correlators: 

n v-A, 2 , 2 2, (0 I J? sl I D s i(p',s))(D sl (p',E) \ J^ A \ B c (p))(B c (p) \ J B ^ | 0) 
[P,P ~ (P' 2 - ™ 2 Dsl )(p 2 - m|J 

+ higher resonances and continuum states , 



A' B ^ D ^(q 2 ) = 


2A B ^ D ^(q 2 ) 
m Bc + m Dsl 


A' B c 
^0 




X B ^ D ^{q 2 ) = 


Af^ Dsl (q 2 ) 
m Bc + m Dal 


a'Bc 
A 2 




Ty B ^ D ^(q 2 ) = 


-2T B ^(q 2 ) , 


T'Bc- 




T' B ^ D ^{q 2 ) = 


_ T Bc ~* Dsl (q 2 ). 
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T]T . PT( 2 l2 2) (0 | J?« | D sl (p',e))(D sl (p',e) \ Jj^_ \ B c (p))(B c (p) | | 0) 
V ^ ) (p' 2 -ml 3l )(p*-ml c ) 

+ higher resonances and continuum states . (8) 

The following matrix elements are defined in the standard way in terms of the leptonic 
decay constants of the D s i and B c mesons as: 

(0 I J u Dsl I D sl (p',e)) = f Dsl m Dsl e» , (0 | J Bc \ B c (p)} = t^^. (9) 

Using Eq. (5), Eq. (6) and Eq. (9) in Eq. (8) and performing summation over the 
polarization of D sl meson we obtain: 

-nV-A/ 2 12 2\ fBj™B c /d s i?™D s i \ ■ a' B c ^>D sl / 2\ a 10 

(m b + m c ) (p /2 - m z Dsi ) (p 2 - m% c ) l 
+ A' B ^ Dsl (q 2 )g^ + A[ B ^ Dsl (q 2 )P^p u + X 2 B ^ Dsl (q 2 )q^Pu\ + excited states, 

T]T-PT ( 2 12 _ fB c m 2 Bc fDsi m Dsi \ T 'B c ^D sl( 2x a //3 

- i T^ Dsl (q 2 )g^ - i T[ B ^ Dsl (q 2 )q^p u ] + excited states. (10) 

To calculate the form factors, A' v , A' , A[, A' 2 , T v , T' and T[, we will choose the 
structures, ie IJiVa pp a p' p , g^, P^Pu, q^Pu, horn nj~ A and e^ a/ 3p a p /l3 , ig^ u and iq^p,, from 
IPj~ PT , respectively. 

On the QCD side, using the operator product expansion (OPE), we can obtain the 
correlation function in quark-gluon language in the deep Euclidean region where p 2 <C 
(mb + m c ) 2 and p <C (m 2 + m 2 ). For this aim, the correlators are written as: 

n y " V, v'\ q 2 ) = * n v - A ^p a p" + n y ~ V + nr- A P M p, + n y -%^, 

^ PT (p 2 ,p' 2 ,q 2 ) = ^ PT e^p a p' P -i^- PT 9, v -i^- PT q^ (11) 

where, each IT, function is defined in terms of the perturbative and non-perturbative parts 
as: 

ITG.V 2 , q 2 ) = Ur( P 2 ,p' 2 , q 2 ) + Ur Per (p 2 , P ' 2 , q 2 ) • (12) 

To obtain the perturbative part of the correlation function, we should study the bare 
loop diagrams in Fig. 1. In calculating the bare loop contributions, we first write the 
double dispersion representation for the coefficients of the corresponding Lorentz structures 
appearing in each correlation function, as: 

nf r = -ttt^ / ds> f ds , Pi{S 2 lf) q2) l2 , + subtraction terms. (13) 
(2tt) 2 J J (s - p 2 )(s' - p' 2 ) 

The spectral densities p? er (s,s',q 2 ) are calculated by the help of the Gutkosky rules, i.e., 
the propagators are replaced by Dirac-delta functions 

1 



p2 _ m 2 



-2in5(p 2 - m 2 ) , (14) 



expressing that all quarks are real. Note that, there are two main vertexes related to 
the bare loop diagrams that describe b — > s transition in Fig. 1, i.e., 7^(1 — 75) 
and a iiV q v {\ + 75). First, we calculate the spectral densities related to 7^(1 — 75) vertex. 
Straightforward calculations end up in the following results: 

Py~ a = 4N C I Q (s, s',q 2 ) {Bi(m b - m c ) - B 2 {m s + m c ) - m c } , 

Po~ A = -2N C I (s, s', q 2 ) {A(m c + m s ) - A'(m b - m c ) - AA±{m b - m c ) 
+2m 2 c (m b — m c — m s ) + m c (2m b m s — u)} , 

p\~ A = 2N C I (s, s\q 2 ){Bi{m b - 3m c ) - B 2 {m c + m s ) + 2A 2 {m b - m c ) 
+2A 3 (m b - m c ) - m c } , 

p\- A = 2N C I (s, s',q 2 ){2A 2 (m b - m c ) - 2A 3 (m b - m c ) - B x {m b + m c ) 

+B 2 (m c + m s ) + m c } . (15) 

Then, the spectral densities related to <J^ u q u (l + 75) vertex are presented as: 

2N C I (s, s', q 2 ) {m c {m b - m s ) + B x [A - m c m s - m\ + m b (m c + m s ) - s] 
B 2 [A' - m c m s - m 2 c + m b (m c + m s ) - s'] }, 



p T v - PT (s,s',q 2 ) = 
+ 



P^ PT (s,s\q 2 ) 



= 2N C I (s, s', q 2 ) {2Ai(2s — u) — A [m b (m c + m s ) — m c {m c + m s ) + s'\ 
+ A' [m b (m c + m s ) — m c {m c + m s ) + s] — 2m c {m b — m c )s' 
— 2m c {m c + m s )s + 2m c {m b + m s )u}, 



p\ PT (s, s', q 2 ) = 2N C I (s, s', q 2 ) {m c (m b + m s ) + B 2 [(m b - m c )(m c + m s ) + s'\ 
+ Bi [A — A' — m b {m c + m s ) + m c {m c + m s ) — s] 

+ (A 2 -A 3 )(2s-u)}, (16) 



where 
I (s,s',q 2 ) = 



1 



A(a, 6, c) = a 2 + 6 2 + c 2 - 2ac - 26c - 2ab, 

A' = (s' + m 2 -m 2 ), 

A = (s + m 2 -m 2 ), 

■u = s + s' — q 2 , 

* 2 = A^V)! 2sA '- Atl l' 

^1 = - „„ ~, ,A (tes'm 2 c - sA' 2 - s'A 2 - t. 2 m 2 + uAA')], 
2\(s,s',q 2 ) 



A 2 = — ——. — -l8ss ,2 m 2 c -2ss'A' 2 -6s' 2 A 2 -2u 2 s'm 2 c 

X 2 (s,s',q 2 ) L 

+ 6s'uAA' -u 2 A' 2 ], 

A 3 = — ■. — Uss'um 2 + 4ss'AA' - 3suA' 2 - 3uA 2 s' - u 3 m 2 + 2h 2 AA'1 . 

A 2 (s, s',q 2 J 

and N c = 3 is the color factor. 

The integration region in Eq. (13) is obtained requiring that the argument of three 
delta vanish, simultaneously. The physical region in the s and s' plane is described by the 
following inequalities: 

2ss' + (s + s' - q 2 )(m 2 - s - m 2 c ) + (m 2 c - m 2 )2s 

A 1 /2(mg,s,m2)A 1 /2( s , s ',g2) ^ 

From this inequality, to use in the lower limit of the integration over s in continuum sub- 
tractions, it is easy to express s in terms of s', i.e., sl is as follow: 

(m 2 c + q 2 -m 2 b - s')(m 2 b s' - q 2 m 2 c ) 

Sl ~ 7 2 2\7 2 A • y l °> 

[mi — q 2 )(m 2 — s'J 

3 Gluon condensate contribution 

In this section, the non-perturbative part contributions to the correlation function are 
discussed. Here, we will follow the same procedure as stated in [38-40,42]. The non- 
perturbative part contains the quark and gluon condensate diagrams. For this aim, we 
consider the condensate terms of dimension 3,4 and 5. It's found that the heavy quark 
condensate contributions are suppressed by inverse of the heavy quark mass and can be 
safely omitted. The light s quark condensate contribution is zero after applying the dou- 
ble Borel transformation with respect to the both variables p 2 and p , because only one 
variable appears in the denominator. Therefore in this case, we consider the two gluon con- 
densate diagrams with mass dimension 4 as non-perturbative corrections. The diagrams for 
contribution of the gluon condensates are depicted in Fig. 2. To obtain the contributions 
of these diagrams, the Fock-Schwinger fixed-point gauge, x^A^ = 0, are used, where A^ is 
the gluon field. In the evaluation of diagrams in Fig. 2, integrals of the following types are 
encountered. 

I (a, 6, c) 
I?(a,b, c) 



(2tt) 4 [jfe2 


— m 2 ] a [(p + k) 2 — m 2 


\ b [(p' + k) 2 - m 2 s } c 


d A k 


k/j, 




(2tt) 4 [k 2 


— m 2 ] a [(p + k) 2 — m 2 


\ b [{p> + k) 2 - ml) c 


d 4 k 


k^k v 




(2tt) 4 [k 2 


— m 2 ] 0, [{p + k) 2 — m 2 


\ b \{p> + k) 2 - m 2 X 



Iuu(a,b,c) = / — — ^— r . (19) 

MV J (2tt) 4 [k 2 - m 2 ] a [{p + k) 2 - m 2 } b \(p' + k) 2 - m 2 ] c 

These integrals can be calculated using the Schwinger representation for the Euclidean 
propagator 

1 1 



(k 2 + m 2 ) n T(n) Jo 



daa n - l e- aik2+m ^ , (20) 

Jo 
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y l P-y 5 )^^V+y 5 ) 
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Y,( 1 -Y 5 )' ,vq v ( 1+ Y5 




y v y 5 





Figure 2: Contribution of gluon condensates for B c — > D sl transition. 

After Borel transformation using: 

B P 2(M 2 )e- ap2 =5(l/M 2 -a) , 

we obtain 



(M 1 2 ) 2 - a - 6 (M 2 2 ) i - a - c W (a + 6 + c -4,l-c-6), 



(21) 



i67r 2 r(a)r(6)r(c) 

4(a,6,c) = A(a,6,c)p^ + I 2 (a,b,c)p ll , 

7^(a, 6, c) = J 6 (a, 6, c)^ + J 3 (a, 6, c)p M ^ + J 4 (a, 6, c)p^ + J 4 (a, 6, c)pjj>„ + I 5 (a, 6, c)p^ . 

(22) 

I in Eq. (22) stands for double Borel transformed form of Eq. (19). In Schwinger repre- 
sentation: 

( _iy+b+c+i 



I k (a,b,c) = i 
I m (a,b,c) = i 
i 6 (a,b,c) = i 



16tt 2 T(a)T(b)T(c) 

i67r 2 r(a)r(6)r(c) 
327r 2 r(a)r(6)r(c) 



(M 1 2 ) 1 - a - 6+fc (M 2 ) 4 ~ a - c - fe W (a + b + c-5,l-c-b) , 
(M*)- a - b - 1+m (M%) 7 - a - c - m U (a + b + c-5,l-c-b) , 
{Mlf- a -\Mlf- a - c U {a + b + c-Q,2-c-b) . (23) 



where k = 1,2, m = 3, 4, 5, M\ and Mf are the Borel parameters in the s and s' channel, 
respectively, and the function U (a, b) is defined as 



where 



U (a,b) 
B 1 



dy(y + M{ + M 2 ) a y b exp 



B-! 



-B - B lV 



1 



M\Ml 
1 

mJm! 

2 

M-M.j 



m\M{ + m 2 b M^ + M 2 Mf(m 2 b + mj - q 2 ) 



[ml + m 2 c ) Ml + Ml (ml + m 2 c ) 



(24) 



Performing the double Borel transformation over the variables p 2 and p' 2 on the physical 
as well as perturbative parts of the correlation functions and equating the coefficients of 
the selected structures from both sides, the sum rules for the form factors A, 



are 



obtained: 



'B C ^D„ 



(m b + m c 



fB c m 2 B J Dsl m Dsl 

CY~ A 



M z M z i 

e i e 2 { — 



4vr 2 



J mj. J St, 



s,s',q 2 )e^e^ 



- %m\m\ (-G 2 

IT 



6 



(25) 



where i = V, 0, 1, 2 and for the form factors Tj 



Bc^D, 



1 , we get 



B C ^D 3 



fB c m 2 B J Dsl m Dal 



M z M z l 

e M i e 2 { — 



An 2 



'°ds> ; 

sl 



pJ- pT (s,s>,q 2 )e^e^ 



r 2 A/r 2 I a s^2\ °jr 



- iMfM 2 (—G 



7T 



T-PT ■ 



6 



(26) 



where j = V, 0, 1. The Sq and s' are the continuum thresholds in B c and -D s i channels, 
respectively and lower bound sl in the integrals is given in Eq. (18). We present the explicit 
expressions of the coefficients C^ A( * T ~ PT ^ correspond to gluon condensates in Appendix-A. 

Now, the ^ Bc ^ D3ll(2) and 7\ Bc ~* Dsll{2) form factors are obtained from the above equa- 
tions replacing the fo sl by decay constant /d s i1(2), and m Dal with m,D al i(2), i- e -> 



A 



'B c ->D al l{2) 



(m b + m c ) 



fBm B J Dsl n2)m Dall{2 ) 
in \ r v ~ A 

2 n /rl / °s ri2 \ °< 



e M i e 



V-A, 



- iM 2 M 2 {^G 



s, s\q 2 )e^e^ 
(27) 
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and 



2 m 2 



jB c m z B J Dsim m Dsim { An 2 Jm% Js L 

qT-PT 



l „2\ „ M'i „ M 2 



- iMfMl^G 2 }^— | . (28) 

Using the straight forward calculations, the form factors of the j\Bc—D sl (2im) anc [ ^ Bc_+Dsl ( 2460 ) 
are found as follows: 

^d s1 (2460) = ( m Bc + m Dsll \ Dsil + / m ge + m Dsl2 \ gsi2 

V m ^ c + m ^ sl / V m ^ c + % sl / 

a b^d si(2 460) = / m^+m^\ ^ + / m Bc + m Dsl \ ^ Dsi2 ^ 

\m Bc + m Dsll J \m Bc +m Dsl2 J 
T B^ Dsl{ 2mo) = ( m Bc +m Dall \ „ Dail ^ + / m Be +m . l2 \ Bc ^ sl2 ^ 

T B C ^D S1 (2460) = T ^ 3l l s ^ +T B C ^ 31 2 C0S ^ (2Q) 

where i' = V,1,2 and / = V, 1. Note that, the Af c ^ Dsl{25m) and y^-*^ 2536 ) form factors 
are obtained from the above equations by replacing the sin9 s — > cos# s and cos# s - * —sin9 s . 



4 Decay widths 

Now, we present the dilepton invariant mass distribution for the B c — > D s \vv and _B C — > 
Dsi/f decays. Using the parameterization of the 5 C — > D s \ transition in terms of form 
factors and also Eq. (3), the dilepton invariant mass distribution of the B c — > D s \vv decay 
can be written as [43]: 



dT (B c -> D sl uu) 3G 2 F m Bc \V tb V;fa 2 \D(x t )\ 2 , 
rfg 2 " 2 8 7r 5 sm 4 w ^ 



SQ!i + -g^Pl 



(30) 



where s = q 2 /m Bc , x t = rn^/m^. The parameters D(x t ), 4> D x , cti and /3i are defined by 



L>(x t ) 



f 2 + x t 3x t — 6 



x t - 1 (x t - l) 2 



In 



(l-r Dsl ) 2 -2 S (l + r Dsl ) + S 2 



(1 - v^T) 2 
(i-v^T) 2 



4rn 



^0 



aB c ^D s1 
^0 



+ 



/\Bc^D sl 



+ 



(i + v^T) 2 

2 V r DA ) 1 + 



+ 



'si 

1 (I -s 



4rD sl (l + v^T) 5 



(31) 
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where td b1 = rn 2 Dsi /m 2 Bc . The differential decay rates for B c — > D sl ll are found to be 



[44,45]: 

dT (B+ - D si m~) G 2 F m 3 Bc \V tb V*\ 2 a 2 , 



dq 2 



2% 5 



-V 0'i 



2m 



1 + — l - )[s a 3 + -^[3 3 )+US 



(32) 



where t = mf /m 2 Bc , v = \Jl — 4m 2 / q 2 and the expressions of 0:3, ^3 and 5 are given as: 



«3 = 



and 
5 = 



(1 " v^) 2 

(i-V^T) 



*^>Q sin 



+ 



nB c ^D s 



y 9 st-V 



+ 



C W A 



(33) 



4r 



s 



r eff a 



Mai _ 2m b C e 7 ff (l + ^¥bZ)n 



+ 



r< /\B c ^D a i 



(1 + v^T) 2 
4r Dsl (l + v^d^; 



z^e// aB c -*D. 
°9 



1 2m fe C 7 e// (l + v/r^DT* 



B c ^D s j 



+ 



r< aBc^Dsi 



r eff A B^ Dsl 2m b C e 7 ff (l + ^b-)Tj 

2m b C 7 e// (l + VT^Ttf 



iB c ^D a 



2 1 r Dil ; 1 + v /5^7 



2 




+ 


^10^1 



B C ^D S 



X 



2m b C e / f (l + ^rb-;)Tl 



B c ^D a 



\C w \ 2 Re(A^ D ^Af^ D ^) \ . (34) 
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20 Dai |A^f - 3(1 - r Dal f\A^ Dsl \ 2 



2(1 + ^)% 
+ -^- [2(1 + r D J - s] \Af^\ 2 + tEnl\ A *^\ 2 

, 0n , (1 - ?d s 



2r 



(35) 



where m?, = m b /mB c 



5 Numerical analysis 

In this section, we present our numerical analysis of the form factors Ai , {i — V, 0, 1, 2) and 
Tj , (j = V, 0, 1). From the sum rules expressions of the form factors, it is clear that the 
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main input parameters entering the expressions are gluon condensates, elements of the CKM 
matrix V t b and V ts , leptonic decay constants f Bc , /d s1 , /d s1 i and /d s1 2, Borel parameters 
Mj 2 and M| as well as the continuum thresholds so and s . We choose the values of the 
condensates (at a fixed renormalization scale of about 1 GeV), leptonic decay constants , 
CKM matrix elements, quark and meson masses as: (^-G 2 ) = 0.012 GeV 4 [46], | V t b \ — 
0.7718:24, I V u |= (40.6 ± 2.7) x lO" 3 [47], Cf f = -0.313, C 9 e// = 4.344, C 10 = -4.669 
[48,49], f Dsl = (225 ± 25) MeV, f Dsll = (240 ± 25) MeV, f Dsl2 = (63 ± 7) MeV [36,50], 
f Bc = (350 ±25) MeV [51-53], m s (l GeV) = (104lgj) MeF, m c = (l. 271^1) CeV, m 6 = 
(4.7±0.07) GeV, m Dsl(2m = (2459.6±0.6) MeV, m Dsl(2536) = (2535.35±0.34±0.5) MeV 
and m Bc = (6.276 ± 0.004) GeV [54]. 

The sum rules for the form factors contain also four auxiliary parameters: Borel mass 
squares M\ and Mf and continuum thresholds s and s' . These are not physical quantities, 
so the the form factors as physical quantities should be independent of them. The parame- 
ters So and s' , which are the continuum thresholds of B c and D s i mesons, respectively, are 
determined from the condition that guarantees the sum rules to practically be stable in the 
allowed regions for M\ and Mf . The values of the continuum thresholds calculated from 
the two-point QCD sum rules are taken to be sq = (45 — 50) GeV 2 and s' = (6 — 8) GeV 2 
[46,55,56]. The working regions for Mf and Mf are determined requiring that not only 
the contributions of the higher states and continuum are small, but the contributions of 
the operators with higher dimensions are also small. Both conditions are satisfied in the 
regions 10 GeV 2 < M 2 < 25 GeV 2 and 4 GeV 2 < Mf < 10 GeV 2 . First, we would like to 
consider the D s \ meson as the pure \cs) state. The values of the form factors at q 2 = are 
presented in Table 2. The sum rules for the form factors are truncated at about 10 GeV 2 , so 



^B c -»D s i(2460)^ 


-0.23 ±0.07 


^B c ->D s i(2536)^ 


-0.22 ±0.06 


^B c ^£> s i(2460) ^ 


0.09 ±0.02 


^B c -+D sl (2536) ^ 


0.07 ±0.02 


jjBc-^D s i(2460) ^ 


0.16 ±0.05 


^S c ^£> 3 i(2536) ^ 


0.17 ±0.05 


jjB c -* D s i(2460)/qn 


-0.26 ±0.08 


^S c ^£> 3 i(2536)/gN 


-0.28 ±0.09 


^B c ^D s i(2460)^Q2^ 


0.12 ±0.03 


rpB c -^D sl (2536) ^ 


0.14 ±0.04 


rpB c ->D al (2460) ^ 


0.11 ±0.03 


rpB c ~* D sl (2536) ^ 


0.14 ±0.04 


rp B c ^D sl (2460) ^ 


-0.14 ±0.04 


rpB c -^D sl (2536) ^ 


-0.16 ±0.05 



Table 2: The value of the form factors of the B c — > D sl (2460) and B c — > D sl (2536) 
transitions at q 2 = 0, Mf = 15 GeV 2 and Mf = 8 GeV 2 , when the D sl meson are 
considered as the pure \cs) state . 

to extend our results to the full physical region, < q 2 < (mg c — m£) gl ) 2 GeV 2 , we look for 
a parameterization of the form factors in such a way that in the region < q 2 < 10 GeV 2 , 
this parameterization coincides with the sum rules predictions. Our numerical calculations 
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show that the sufficient parameterization of the form factors with respect to q 2 is as follows: 

/.( g 2 ) = + b — — . (36) 

(I — q , \ (I — -t-V 

v mi., i v .. / 

/at /it 

The values of the parameters a,b and to/^ are given in the Table 3. To calculate the 





a b rrifit 




a b rrifit 


A B C ->D S1 (2460) ( 


-0.13 -0.10 5.30 


^B c -D s i(2536)/ 2 \ 


-0.12 -0.10 5.22 


A B c ^ sl (2460) (g2) 


0.05 0.04 5.98 


^B c -D s i (2536)^2) 


0.04 0.03 5.99 


^B c -D s i (2460)^2) 


0.09 0.07 5.95 


^B c -D s i(2536)^2) 


0.09 0.08 5.98 


^c-D sl (2460) (?2) 


-0.16 -0.10 5.15 


^B c -D s i(2536), 2 \ 


-0.17 -0.11 5.17 


,pB c ->D s i(2460)^2) 


0.08 0.04 5.22 


y£ c —D s i(2536)^ 


0.09 0.05 5.05 


,pB c ->D s i(2460)^2) 


0.10 0.01 5.85 


rpB c ^D sl (2536) 


0.11 0.03 5.96 


yS c ->D s i(2460)^ 2 ^ 


-0.09 -0.05 5.22 


rpB c ->D sl (2536) ^2^ 


-0.08 -0.07 5.14 



Table 3: Parameters appearing in the fit function for the form factors of the B c — > _D s i(2460) 
and B c -> £> sl (2536) transitions at M 2 = 15 GeV 2 and Mf = 8 GeV 2 . 

branching ratios of the B c — > /} sl (2460[2536])Z + Z~/W decays, we integrate Eqs. (30, 32) 
over q 2 in the whole physical region and use the total mean life time tb c = (0.46 ± 0.07) ps 
[54] . Our numerical analysis shows that the contribution of the non-perturbative part (the 
gluon condensate diagrams ) is about 12% of the total and the main contribution comes 
from the perturbative part of the form factors. The values for the branching ratio of these 
decays are obtained as presented in Table 4, when only the short distance (SD) effects are 
considered. It should be noted that, the long distance (LD) effects for the charged lepton 



MODS 


BR 


MODS 


BR 


B c -> D sl (24Q0)uu 


(3.26 ± 1.10) x 10~ 7 


B c -> D sl (2536)uP 


(2.76 ±0.88) x 10~ 7 


B c -> D s i(2460)e + e" 


(5.40 ± 1.70) x 10~ 6 


B c -> D sl (2536)e+e- 


(2.91 ±0.93) x 10~ 6 


B c -> D sl (2460)/i+/x- 


(2.27 ±0.95) x 10~ 6 


B c -> D sl {2536)fi + fi- 


(1.96 ±0.63) x 10- 6 


B c ^ D sl (2A60)r + r- 


(1.42 ±0.45) x 10~ 8 


B c -> D sl (2536)r+r- 


(0.68 ±0.21) x 10~ 8 



Table 4: The branching ratios of the semileptonic B c — > D sl (2460)/ + / jvv and B c 
D sl (253G)l + l~ jvv decays with SD effects. 
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modes are not included in the values of Table 4. With the LD effects, we introduce some 
cuts close to q 2 = and around the resonances of J ftp and ip' and study the three regions 
as follows: 

i ■ 4vLn < \[q 2 < M m - 0.20, 

II : Mj/^ + 0.04 < \fq 2 < M^i — 0.10, 
///: ilfy +0.02 < yg2 < m Bc -m Dal . (37) 



where yq^in = ^ m i- I* 1 Table 5, we present the branching ratios in terms of the regions 
shown in Eq. (37). The errors are estimated by the variation of the Borel parameters M\ 



MODS 


7 


// 


III 


BR[B C -> D s i(2460)e+e-) 


(4.40 ± 1.35) x 10^ 6 


(1.62 ±0.52) x 10~ 7 


(1.01 ±0.35) x 10" 7 


BR(B C -> D s i(2536)e+e-) 


(2.65 ±0.82) x 10" 6 


(0.90 ±0.28) x 10~ 7 


(0.81 ±0.25) x 10- 7 


BR(B C -> D sl (2460)^ + Ai") 


(1.76 ±0.58) x 10- 6 


(1.61 ±0.53) x 10~ 7 


(1.01 ±0.35) x 10- 7 


Bi2(B c -> D sl (2536)/i+/i^) 


(1.00 ±0.31) x 10- 6 


(0.96 ±0.29) x 10- 7 


(0.68 ±0.21) x 10~ 7 


Bi2(B c -> D sl (2460)r+r-) 


undefined 


(5.20 ± 1.61) x 10" 9 


(8.12 ±2.75) x 10" 9 


£i2(£ c -> D sl (2536)r+r-) 


undefined 


(4.16 ± 1.29) x 10~ 9 


(6.56 ±2.16) x 10~ 9 



Table 5: The branching ratios of the semileptonic B c — > D s i(2460)/ + / jvv and B c — > 
D s i{2536)l + l^ /uu decays including LD effects. 

and M|, the variation of the continuum thresholds s and s' Q , the variation of b and c quark 
masses and leptonic decay constants fs c and fo sl - 

Now, we would like to analyze the form factors obtained when we considered the D s i 
meson as a mixture of two \ 3 Pi) and I 1 .Pi) states (see Eq. 29). The transition form factors 
of B c -> D sl (24Q0[2536])l + l- /vis at q 2 = in the interval -180° < 9 S < 180° are shown in 
Figs. 3 and 4. From these figures, we see that all form factors have the following common 
behaviors: 1) they have extrema at the same mixing angles and 2) they come across at 
two points. The dependence of the form factors B c — > D sl (2460) and B c — * D sl (2536) 
transitions on the mixing angle, S , and the transferred momentum square, q 2 , are plotted 
in Figs. 5, 6 in the regions < q 2 < {m Bc — m Dal ) 2 GeV 2 and —180° < 6 S < 180° for q 2 
and mixing angle, respectively. Using Eqs. 30 and 32 we analyze the decay widths and the 
branching ratios related to considered decays. For this aim, we denote the variations of the 
decay widths with respect to q 2 and 9 S in the regions 4mf < q 2 < {ms c — rriD sl ) 2 GeV 2 , 
(I = r, fi, e) and —180° < 9 S < 180° and branching ratios only in terms of mixing angle 9 S in 
Figs. 7-10. The results for electron and muon are approximately the same, so we consider 
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I — /jl,t. In Fig. 11, as an example, we only depict the variation of the branching ratio of 
B c — > D sl (2460)uu decay in terms of the mixing angle. The figures 9 and 10 depict a regular 
variation of the branching ratios for case with respect to the mixing angle, while we 
see an irregular variation of the branching ratio of the B c — > D s i{2AQQ)uu transition with 
respect to the 9 S . 

In summary, We analyzed the semileptonic B c — > D sl (2460[2536])/ + / _ , / = e,/i,r and 
B c — > D sl (2460 [2536] )vv decays in the framework of the three-point QCD sum rules. First, 
we assumed the -D sl (2460) and _D sl (2536) axial vector mesons as the pure \cs) states. In this 
case, the related form factors were computed. The branching ratios of these decays were 
also estimated with both the short distance (SD) and long distance (LD) effects, for the 
charged lepton modes. Second, D s \ (2460 [2536]) mesons were considered as a combinations 
of two states \ 3 Pi) = |-D s il) and I 1 Pi) = \D s i2) with different masses and decay constants. 
We evaluated the transitions form factors and the decay widths of these decays with respect 
to the mixing angle 9 S and the transferred momentum square q 2 . The dependence of the 
branching ratios on 9 S was also presented. Detection of these channels and their comparison 
with the phenomenological models like QCD sum rules could give useful information about 
the structure of the D s i meson and the mixing angle 9 S . 
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Appendix-A 

In this appendix, the explicit expressions of the coefficients of the gluon condensate entering 
the sum rules of the form factors X.* C ^ D *\ (i = V,0, 1,2) and rj Bc ^ Dsl , (j = V,0, 1) are 
given. 



C 



V-A 
V 



r<v-A 



-10/1(3,2 
+10/1(3,2 

+60/ 2 (l,4 
+60/ (l,4 
+10/i(3,2 

-20 /J ' 11 (3 
-10/o(3,2 

-60/i(l,4 
+10/o(3,2 
-20/i(2,2 
-10/o(3,l 
-20 /J ' 11 (3 
-20/l 0,1] (3 
+ 10/f' 2] (3 
+30/o(2,2 
+20/{°' 1] (3 
+20/ 2 (2,2 
-20/o(l,2 
+20/ 2 °' 1] (2 
-10 / 2 0,2] (3 



2)m b 3 m c 2 + 10 /i(3, 2, 2)m b 2 m c 3 + 10 / 2 (3, 2, 2)m b 2 m c 3 + 10 / (3, 2, 2)m b 2 m c 3 
2)m b m L A + 10 / 0,11 (3, 2, 2)m b 2 m c - 30 /i(3, 2, l)m b 2 m c + 60 /i(l, 4, l)m 6 2 m c 
l)m b 2 m c - 20 / 2 (3, 2, l)m b 2 m c + 10 / 2 °' 1] (3, 2, 2)m b 2 m c - 20 / (3, 2, l)m b 2 m c 
l)m b 2 m c + 10/f' 1] (3, 2, 2)m b 2 m c + 20 /i(2, 2, 2)m fe m c 2 + 10 / 2 (3, 2, l)m fe m c 2 

1) m 6 m c 2 + 40 / 2 (2, 3, l)m b m 2 - 10 / (3, 2, l)m b m c 2 + 20 h{2, 3, l)m 6 m c 2 
2, 2)m 6 m c 2 + 30 h (4, 1, l)m 6 m c 2 - 10 / 2 (3, 2, 2)m c 5 - 10 A(3, 2, 2) 

2) m c 5 + 20 A(3, 2, l)m b 3 + 10 A(2, 2, 2)m b 3 - 20 A(2, 3, l)m b 3 
l)m b 3 - 10 /j 0,11 (3, 2, 2)m b 3 - 30 / 2 (4, 1, l)m c 3 + 20 / 2 0,11 (3, 2, 2)m c : 

1) m c 3 - 10 / 2 (3, 1, 2)m c 3 - 20 / (2, 2, 2)m c 3 - 20 / 2 (2, 2, 2)m c 3 

2) m c 3 - 30 / (4, 1, l)m c 3 - 30 h (4, 1, l)m c 3 + 20 /f' 1] (3, 2, 2)m c 3 
2)m c 3 + 20 / 0,11 (3, 2, 2)m c 3 - 50 /i(2, 2, l)m 6 + 20 /{°' 1] (2, 3, l)m b 
2, l)m b + 20 A(l, 2, 2)m 6 + 60 / (1, 3, l)m b - 20 / 2 (2, 2, l)m b 
1, 2)m fe - 20 / (2, 2, l)m 6 + 30 A(2, 1, 2)m 6 + 100 / 2 (1, 3, l)m 6 



2, 2)m fe - 20 /j 0,1| (2, 2, 2)m fe + 40 / 0,1| (2, 3, l)m b + 20 h(l, 3, l)m 6 
l)m c + 30 / 2 °' 1] (3, 1, 2)m c + 20 / 2 °' 1] (3, 2, l)m c + 10 / °' 1] (3, 2, l)m c 
2, l)m c - 10 / °' 21 (3, 2, 2)m c + 20 /1°' 1] (3, 1, 2)m c + 20 /{°' 1] (2, 2, 2)m c 

1) m c - 30 / 2 (2, 1, 2)m c + 10 / (3, 1, l)m c + 20 / °' 1] (2, 2, 2)m c 

2) m c - 20 / 2 (1, 2, 2)m c + 30 / 0,11 (3, 1, 2)m c - 10 /i(3, 1, l)m c 
2, 2)m c - 10 / 2 (3, 1, l)m c - 20 A(2, 1, 2)m c - 30 / (2, 1, 2)m c 

2, 2)m c + 20 /i(2, 2, l)m c - 20 /i(l, 2, 2)m c - 10 /{ 0,2] (3, 2, 2)m c 



-20 / 6 (3, 2, 2)m c 5 - 40 / 6 (3, 2, l)m c 3 - 20 / 6 (3, 1, 2)m c 3 + 40 4°' 6] (3, 2, 2)m c 3 
-40 / 6 (2, 2, 2)m c 3 - 60 / 6 (4, 1, l)m c 3 + 5 / (3, 1, l)m c 3 - 30 / °' 1] (1, 4, l)m fe 3 
+20 / 6 (2, 2, 2)m b 3 + 5 / (2, 2, l)m b 3 - 120 / 6 (1, 4, l)m b 3 + 40 / 6 (2, 3, l)m 6 3 
+ 15 / °' 1] (3, 2, l)m fe 3 - 20 / (1, 3, l)m 6 3 + 10 / °' 1] (2, 3, l)m fe 3 - 5 / °' 2] (3, 2, 2)m b 
+10 /I°' 1] f2. 2, 2)m 6 3 - 5 / (1, 2, 2)m 6 3 - 20 4 * 1] (3, 2, 2)m b 3 + 20 4°' 1] (3, 1, 2)m c 

)m c + 40 /f' 11 (2, 2, 2)m c - 20 4°' 2] (3, 2, 2)m c + 20 / 6 (3, 2, l)m 6 3 
+5 / ' 1] (3, 1, l)m c + 5 / (1, 1, 2)m c + 20 / 6 (2, 1, 2)m c + 40 / 6 (3, 1, l)m c 
+5 / (1, 2, l)m c - 5 / (2, 1, l)m c - 10 / °' 2] (2, 3, l)m 6 - 10 / ' 2] (3, 1, 2)m 6 
-10 / 0,11 (1, 3, l)m b - 15 / (1, 2, l)m b - 40 / 6 (2, 2, l)m b + 15 / °' 1] (2, 2, l)m b 
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+20 /{ 0>1I (1, 2, 2)m b - 40 7 6 (1, 3, l)m b - 40 7 6 (3, 1, l)m b + 40 2, l)m c 

-20 /| 0,11 (2, 2, 2)m fe + 20 /|°* 2l (3, 2, 2)m fe - 40 /| 0,11 (3, 1, 2)m 5 - 15 J (l, 1, 2)m 5 
-15 4°* 21 (2, 2, 2)m fc + 15 7 (2, 1, l)m 6 - 20 7 6 (2, 1, 2)m b + 25 j£' 1] (2, 1, 2)m b 
+ 10 4°' 1] (3, 1, l)m b - 15 7 °' 2] (3, 2, l)m b - 20 7 6 (1, 2, 2)m b - 40 >1] (2, 3, l)m b 
-20 (3, 2, l)m b - 5 7 (3, 2, 2)m c 6 m b + 10 7 (2, 2, 2)m c 3 m b 2 + 20 7 6 (3, 2, 2)m 3 m b 2 
-10 7 (2, 3, l)m c 4 m fe + 15 4°' 1] (3, 2, 2)m c 4 m b + 20 7 6 (3, 2, 2)m c 4 m 5 - 15 7 (2, 2, 2)m c 4 m fe 
-5 7 (3, 2, l)m c 4 m 6 - 15 7 (4, 1, l)m c 4 m 6 - 5 7 (3, 2, 2)m c 3 m 6 4 + 5 7 (3, 2, 2)m c 4 m 6 3 
+5 7 (3, 2, 2)m c 5 m b 2 - 30 7 (1, 4, l)m c m fe 4 - 5 4°' 1] (3, 2, 2)m c m b 4 + 10 7 (3, 2, l)m c m fe 4 
+ 10 7 (2, 3, l)m c 2 m b 3 - 10 7 (3, 2, l)m c 2 m b 3 - 20 7 6 (3, 2, 2)m 2 m b 3 + 30 7 (1, 4, l)m 2 m b 3 
-10 4 0,1] (3, 2, 2)m 3 m b 2 + 5 7 (3, 2, l)m 3 m b 2 + 15 7 (4, 1, l)m 3 m b 2 + 20 7 6 (2, 2, 2)m 2 m b 
+ 10 4 ' 1] (3, 1, 2)m c 2 m 6 - 40 (3, 2, 2)m 2 m b - 15 4°' 2] (3, 2, 2)m 2 m b 
+10 7 (1, 3, l)m c 2 m fe + 20 4° ,1] (3, 2, l)m c 2 m fe - 20 7 (1, 2, 2)m c 2 m 6 - 15 7 (2, 1, 2)m c 2 m fc 
-40 7 6 (2, 3, l)m 2 m b + 15 4 ' 1] (4, 1, l)m 2 m b + 30 j{ 0,11 (2, 2, 2)m c 2 m 6 + 60 7 6 (4, 1, l)m 2 m b 
-10 7 (3, 1, l)m 2 m b + 20 7 6 (3, 1, 2)m c 2 m 6 + 15 7 (2, 2, l)m c 2 m b + 20 4° ' 1] (2, 3, l)m 2 m b 
+40 7 6 (3, 2, l)m c 2 m 6 - 10 4°' 1] (2, 2, 2)m c m fe 2 - 20 7 6 (3, 2, l)m c m b 2 + 20 4°' I] (3, 2, 2)m c m fe 2 
+15 7 (2, 1, 2)m c m b 2 + 5 7 (3, 1, l)m c m b 2 - 20 4° ,1] (3, 1, 2)m c m 6 2 - 20 7 6 (2, 2, 2)m c m fe 2 
-30 7 (1, 3, l)m c m b 2 + 120 7 6 (1, 4, l)m c m b 2 + 10 7 (1, 2, 2)m c m b 2 - 5 4° ,1] (3, 2, l)m c m b 2 
-10 7 (2, 2, l)m c m fe 2 + 5 4°' 2] (3, 2, 2)m c m b 2 

-40 4 0,11 (2, 3, l)m b + 20 /|°' 2] (3, 2, 2)m fe - 40 7 3 (2, 2, l)m 6 - 20 7\(1, 2, 2)m b 

-20 h{2, 1, 2)m fe - 40 7 3 °' 1] (2, 3, l)m b - 20 I^ 1] (3, 2, l)m b - 20 7 3 (2, 1, 2)m b 

-20 7 3 0,11 (2, 2, 2)m fe - 20 7 3 (1, 2, 2)m b - 20 7 4 (1, 2, 2)m 6 - 10 If ' 1] (2, 3, l)m b 

+5 7{°' 2] (3, 2, 2)m 6 - 5 7 2 (3, 2, 2)m 3 5 - 5 7 (3, 2, 2)m c 5 - 20 7 3 (3, 2, 2)m c 5 

-15 7\(3, 2, 2)m c 5 - 45 4(3, 2, l)m 3 - 20 7 4 (3, 1, 2)m 3 - 20 7 2 (3, 2, l)m 3 

-15 7 2 (4, 1, l)m c 3 - 25 7 (3, 2, l)m c 3 - 40 7 4 (3, 2, l)m c 3 - 10 7 2 (2, 2, 2)m c 3 

-45 A (4, 1, l)m c 3 - 20 7 4 (3, 2, 2)m c 5 - 5 7 (3, 1, 2)m c 3 - 40 7 3 (2, 2, 2)m c 3 

+10 4°' 1] (3, 2, 2)m c 3 - 30 7 4 (2, 2, 2)m c 3 + 40 7 3 0,1] (3, 2, 2)m c 3 - 60 7 3 (4, 1, l)m c 3 

-15 7 (4, 1, l)m c 3 - 5 7 2 (3, 1, 2)m c 3 - 10 7 (2, 2, 2)m c 3 - 20 7 4 (3, 1, 2)m c 3 

-40 7 3 (3, 2, l)m c 3 + 30 7{ 0,1] (3, 2, 2)m c 3 + 40 4 0,11 (3, 2, 2)m c 3 - 60 7 4 (4, 1, l)m c 3 

-20 7 3 (3, 2, 2)m c 2 m fe 3 + 20 7 4 (2, 2, 2)m 2 m b - 20 7 (2, 3, l)m 2 m b + 40 7 4 (3, 2, l)m 2 m b 

-40 4°' 11 (3, 2, 2)m 2 m b - 20 7 2 (2, 3, l)m 2 m b + 60 7 4 (4, 1, l)m 2 m b + 60 7 3 (4, 1, l)m 2 m b 

+20 7 3 (3, 1, 2)m 2 m b + 20 7 3 (2, 2, 2)m c 2 m fe + 5 A(3, 2, 2)m c 4 m fe + 20 7 3 (3, 2, 2)m c 4 m fe 

+20 7 4 (3, 2, 2)m c 4 m fe + 5 7 (3, 2, 2)m 3 m b 2 + 20 7 4 (3, 2, 2)m 3 m b 2 + 5 7 2 (3, 2, 2)m c 3 m fe 2 

+15 7i(3, 2, 2)m c 3 m fe 2 + 20 7 3 (3, 2, 2)m 3 m 2 - 20 7 4 (3, 2, 2)m c 2 m fe 3 - 5 A(3, 2, 2)m c 2 m 6 3 

+ 15 7 2 (3, 2, l)m 2 m b + 15 7 4 (4, 1, l)m 2 m b + 5 7 (3, 2, l)m 2 m b + 10 7 4 (3, 1, 2)m 2 m b 

-50 A(2, 3, l)m c 2 m fc - 10 If' 11 (3, 2, 2)m 2 m b + 35 7i(3, 2, l)m 2 m b + 20 7 4 (3, 1, 2)m 3 2 m b 
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-40 /[ ' 1] (3, 2, 2)m c 2 m b - 40 4(2, 3, l)m c 2 m b - 20 4(3, 1, 2)m c 3 + 10 7 2 °' 1] (3, 2, 2)m c 3 

+40 4(2, 3, l)m fe 3 + 20 4(2, 2, 2)m 6 3 - 5 (3, 2, 2)m 6 3 + 40 4(2, 3, l)m b 3 

-20 /f 1] (3, 2, 2)m 6 3 - 120 4(1, 4, l)m 6 3 + 20 4(3, 2, l)m 5 3 - 120 4(1, 4, l)m b 3 

+20 4(3, 2, l)m 6 3 + 10 4(2, 3, l)m 6 3 - 30 4(1, 4, l)m 6 3 - 40 4(2, 2, 2)m c 3 

+40 4(3, 2, l)m c 2 m b - 40 4(2, 3, l)m c 2 m b + 20 (3, 2, 2)m c m b 2 + 5 4°' 1] (3, 2, 2)m c m fe 2 

-20 4(2,2,2)m c m 6 2 - 30.fi (3, 2, l)m c m b 2 + 90 4(1, 4, l)m c m 6 2 + 120 4(1, 4, l)m 3 m fe 2 

+20 J 3 °' 1] (3, 2, 2)m c m b 2 + 120 4(1, 4, l)m c m b 2 + 10 4(3, 2, l)m fe 3 - 20 J 3 °' 1] (3, 2, 2)m b 3 

+40 4(3, 1, l)m c - 5 4(2, 2, l)m c + 10 4°* 1] (2, 2, 2)m c + 20 4(2, 1, 2)m c + 40 7 3 °' 1] (3, 2, l)m c 

+20 /f 1] (3, 1, 2)m c + 10 4°' 1] (2, 2, 2)m c + 20 J 3 (2, 2, l)m c - 20 /f 2] (3, 2, 2)m c 

+40 J 4 °' 1] (2, 2, 2)m c + 5 J 2 (3, 1, l)m c + 20 J 4 (2, 2, 2)m 6 3 + 30 J (l, 4, l)m c m fe 2 

+5 7i(3, 1, 2)m c m b 2 - 15 7 (3, 2, l)m c m b 2 + 5 4°' 1] (3, 2, 2)m c m fe 2 - 10 7 2 (3, 2, l)m c m b 2 

+30 7 2 (1, 4, l)m c m 6 2 - 20 7 4 (3 , 2, l)m c m fe 2 + 15 7{°' 1] (3, 2, 2)m c m 6 2 - 10 h(2, 2, 2)m c m b 2 

-20 7 3 (3, 2, l)m c m fe 2 - 20 7 3 (2, 2, 2)m c m b 2 - 15 j{°' 2] (3 ? 2 , 2)m c - 10 h(2, 1, 2)m c 

-5 7 2 °' 2] (3, 2, 2)m c + 5 h(2, 2, l)m c + 40 if ,1] (3, 2, l)m c + 10 (3, 2, l)m c 

+15 4°' 1] (3, 1, 2)m c + 20 7 3 °' 1] (3, 1, 2)m c + 35 f} 0,1] (3, 2, l)m c + 30 f} 0,1] (3, 1, 2)m c 

-5 4° ' 2] (3, 2, 2)m c + 40 7 3 °' 1] (2, 2, 2)m c + 20 7 3 (2, 1, 2)m c - 15 7 (2, 1, 2)m c 

-40 7 3 0,11 (3, 1, 2)m b - 40 7 4 (3, 1, l)m 6 - 10 7|°' 1] (3, 2, l)m 6 - 40 7 4 (2, 2, l)m 6 

+20 7 3 °' 2] (3, 2, 2)m 6 - 40 7 3 (1, 3, l)m b - 40 7f 1] (3, 1, 2)m b + 10 7^1, 3, l)m b 

-10 7 (2, 2, l)m b - 20 7f 1] (2, 2, 2)m b - 10 7 2 (2, 2, l)m b - 10 7f' 1] (3, 1, 2)m b 

+10 7 (1, 3, l)m 6 - 20 7f 1] (3, 2, l)m 6 - 20 7 3 °' 2l (3, 2, 2)m c - 10 7 (3, 1, l)m c 

+40 7 4 (3, 1, l)m c - 15 7 2 (2, 1, 2)m c + 20 7 4 (2, 2, l)m c + 15 7i(3, 1, l)m c + 15 7 2 0,11 (3, 1, 2)m c 

+15 If ' 1] (3, 2, l)m c - 40 7 3 (3, 1, l)m 6 - 20 7 4 (2, 1, 2)m 6 + 10 7 2 (1, 3, l)m b - 30 h(2, 2, l)m b 

-40 7 4 (1, 3, l)m b + 30 7{°' 1] (2, 2, 2)m c 

C 2 V '- A = 15 7 2 (4, 1, l)m 2 m b - 40 7 3 ' 1] (3, 2, 2)m c 2 m b - 40 7 4 (3, 2, l)m c 2 m fe - 10 7 2 (2, 3, l)m 3 2 m fe 
+20 7 3 (3, 1, 2)m 2 m b + 60 7 3 (4, 1, l)m 2 m b - 20 7 4 (2, 2, 2)m 2 m b + 40 7 3 (3, 2, l)m 2 m b 
+40 7i°' 1] (3, 2, 2)m c 2 m fe - 60 7 4 (4, 1, l)m 2 m b + 40 7 4 (2, 3, l)m 2 m b + 20 7 3 (2, 2, 2)m 2 m b 
-5 7 (3, 2, l)m 2 m b + 15 7i(3, 2, l)m 2 m b - 20 A(2, 3, l)m 2 m b - 40 7 3 (2, 3, l)m 2 m b 
-20 7 4 (3, 1, 2)m c 2 m fe + 10 7 2 0,1] (3, 2, 2)m c 3 + 60 7 4 (4, 1, l)m c 3 - 20 7 3 (3, 1, 2)m c 3 
+40 7 4 (2, 2, 2)m c 3 - 15 7\(4, 1, l)m c 3 - 40 7 3 (3, 2, l)m c 3 - 40 7 3 (2, 2, 2)m c 3 
-15 7 2 (4, 1, l)m c 3 - 5 7 2 (3, 2, l)m c 3 + 10 7{°' 1] (3, 2, 2)m c 3 + 5 7 (3, 1, 2)m c 3 
-10 4 0,11 (3, 2, 2)m c 3 - 10 7i(2, 2, 2)m c 3 + 10 7 (2, 2, 2)m c 3 - 10 7 2 (3, 1, 2)m c 3 
-5 4(3, 1, 2)m c 3 + 15 7 (4, 1, l)m c 3 - 20 4(3, 2, l)m c 3 + 20 4 ' 1] (3, 2, 2)m 6 3 
-30 4(1, 4, l)m 6 3 - 20 7 3 °' 1] (3, 2, 2)m b 3 + 40 4(2, 3, l)m fe 3 + 5 4(3, 2, l)m 2 m b 
-20 4°' 1] (3, 2, 2)m c m fe 2 + 5 4(3, 1, 2)m c m b 2 - 20 4(3, 2, l)m c m fe 2 + 20 4(2, 2, 2)m c m fe 2 
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-5 (3, 2, 2)m c m b 2 - 10 7 2 (3, 2, l)m c m 6 2 + 15 7 (3, 2, l)m c m fe 2 + 20 7 3 °' 1] (3, 2, 2)m c m b 2 

-10 7 2 (2, 2, 2)m c m fe 2 - 30 7 (1, 4, l)m c m 6 2 + 120 7 3 (1, 4, l)m c m fe 2 + 5 7{°' 1] (3, 2, 2)m c m fe 2 

+30 h{l, 4, l)m c m b 2 - 120 7 4 (1, 4, l)m c m b 2 + 5 7 2 °' 1] (3, 2, 2)m c m 5 2 - 20 7 3 (2, 2 , 2)m c m fe 2 

+20 7 4 (3, 2, l)m c m 6 2 - 10 A(3, 2, l)m c m 6 2 + 30 7 2 (1, 4, l)m c m b 2 + 20 7 3 (2, 2, 2)m 6 3 

-20 7 4 (2, 2, 2)m 6 3 - 5 7 2 ' 1] (3, 2, 2)m fe 3 - 40 7 4 (2, 3, l)m b 3 + 120 7 4 (1, 4, l)m b 3 

-20 7 4 (3, 2, l)m 6 3 + 10 7 2 (2, 3, l)m b 3 + 10 7 2 (3, 2, l)m 6 3 - 120 7 3 (1, 4, l)m fe 3 

+20 7 3 (3, 2, l)m fe 3 + 20 7 4 °' 2] (3, 2, 2)m c + 10 7{ 0,1] (3, 2, l)m c - 40 7 4 °' 1] (2, 2, 2)m c 

+5 7 (2, 2, l)m c + 15 If' 11 (3, 1, 2)m c - 5 Cl 2 (3, 2, 2)m c - 15 (3, 2, l)m c 

+5 7i 0,2l (3, 2, 2)m c - 40 7 4 (3, 1, l)m c + 15 7 (2, 1, 2)m c + 5 7 2 (3, 1, l)m c 

+40 7 3 (3, 1, l)m c + 5 7^3, 1, l)m c - 20 7i°' 1] (3, 1, 2)m c + 10 7 (3, 1, l)m c 

+20 7 3 (2, 2, l)m c + 15 7 2 °' 1] (3, 2, l)m c - 20 7 4 (2, 2, l)m c + 20 7 2 (2, 1, 2)m c 

+10 i|°' 1] (2, 2, 2)m c + 20 7 3 °' 1] (3, 1, 2)m c - 20 7 4 (2, 1, 2)m c + 20 7 3 (2, 1, 2)m c 

-40 7f 11 (3, 2, l)m c + 40 7 3 0,11 (3, 2, l)m c - 5 if ' 2] (3, 2, 2)m c + 40 7 3 ' 1] (2, 2, 2)m c 

+5 7 2 (2,2,l)m c + 10 7 2 0,11 (2,2,2)m c - 15 7!(2,l,2)m c - 15 4 0,1] (3, l,2)m c 

-10 7i 0,1] (2, 2, 2)m c - 20 7 3 0,2] (3, 2, 2)m c + 10 7^1, 3, l)m b + 20 7 4 (1, 2, 2)m b 

-40 7 3 0,1] (3, 1, 2)m b + 40 7f 1] (3, 1, 2)m b - 20 7 2 (2, 1, 2)m 6 + 20 7f 1] (3, 2, l)m b 

-10 7 2 0,11 (3, 2, l)m 6 - 20 7 3 (1, 2, 2)m 6 - 10 7 2 °' 1] (2, 3, l)m b - 40 7 3 (1, 3, l)m b 

+5 7 2 °' 2] (3, 2, 2)m 5 + 10 7 (2 , 2, l)m 6 - 20 7 3 (2, 1, 2)m b - 10 7x(2, 2, l)m 6 

-10 7 2 (2, 2, l)m 6 - 10 7 (1, 3, l)m b - 20 7 4 °' 2] (3, 2, 2)m 6 - 40 7 3 °' 1] (2, 3, l)m b 

-20 7 2 (1, 2, 2)m 6 + 20 7 4 (2, 1, 2)m b - 40 7 3 (3, 1, l)m 6 - 10 7 2 (1, 3, l)m b 

-10 7 2 0,1] (3, 1, 2)m b - 20 7 3 0,1] (2, 2, 2)m 6 + 40 7 4 (1, 3, l)m b + 40 7 4 (3, 1, l)m b 

+20 7 3 °' 2] (3, 2, 2)m fe - 40 7 3 (2, 2, l)m 6 + 20 7 4 °' 1] (2, 2, 2)m fe + 40 7 4 °' 1] (2, 3, l)m b 

-20 7 3 °' 1] (3, 2, l)m 6 + 40 7 4 (2, 2, l)m 6 + 25 7 (3 , 2, l)m c 3 - 60 7 3 (4, 1, l)m c 3 

+10 7 2 (3, 1, 2)m c 2 m b + 20 7 4 (3, 2, 2)m c 5 + 5 7 (3, 2, 2)m c 5 - 5 7 2 (3, 2, 2)m c 5 

-20 7 3 (3,2,2)m c 5 - 5/i(3, 2, 2)m c 5 - 40 If 1] (3, 2, 2)m c 3 - 10 7 2 (2, 2, 2)m c 3 

+40 7 4 (3, 2, l)m c 3 + 40 7 3 0,1] (3, 2, 2)m c 3 + 20 7 4 (3, 1, 2)m c 3 + 5 7 2 (3, 2, 2)m c A m b 

-20 7 4 (3, 2, 2)m c A m b + 20 7 3 (3, 2, 2)m c A m b + 5 7i(3, 2, 2)m c 3 m b 2 - 5 7 (3, 2, 2)m c 3 m b 2 

+20 7 3 (3, 2, 2)m c 3 m b 2 - 20 7 4 (3, 2, 2)m 3 m b 2 + 5 7 2 (3, 2, 2)m c 3 m b 2 - 5 7 2 (3, 2, 2)m 2 m b 3 

+20 7 4 (3, 2, 2)m c 2 m fe 3 - 20 7 3 (3, 2, 2)m c 2 m b 3 - 10 7 2 0,11 (3, 2, 2)m 2 m b + 20 7 (2, 3, l)m 2 m b 

d~ PT = 10 A (3, 2, 2)m b 4 m c 2 - 20 7\(1, 1, 2) + 10 4°' 1] (2, 2, 1) - 10 7 2 °' 2] (3, 2, 1) 

-10 7 (3, 2, 2)m fe 3 m c 3 - 10 A (3, 2, 2)m b 3 m c 3 - 10 ^(1, 2, 1) - 10 7 (1, 1, 2) 
-60 4, l)m fe 3 m c - 10 4°' 1] (3, 2, 2)m b 3 m c - 60 7 (1, 4, l)m fe 3 m c + 20 7 2 (3, 2, l)m fe 3 m c 
+20 7i(3, 2, l)m 6 3 m c - 10 7 2 0,1] (3, 2, 2)m fe 3 m c - 10 7{°' 1] (3, 2, 2)m 6 3 m c - 60 7 2 (1, 4, l)m 6 3 m, 
+20 7 2 (2, 3, l)m b 3 m c - 10 7 2 (3, 2, 2)m 6 3 m c 3 + 10 7 2 (3, 2, 2)m b m c 5 + 10 7i(3, 2, 2)m fe m c 5 
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+10 J„(3, 2, 2)m b m c 5 + 20 J (3, 2, l)m 6 3 m c - 20 h(2, 3, l)m b 3 m c + 10 A(3, 1, 2)m b 3 m c 

+20 7 (2, 3, l)m b 3 m c - 10 h(3, 2, 2)m fe 2 m c 4 + 10 (3, 2, 2)m 6 4 - 20 h(3, 2, l)m fe 4 

+60/i(l,4,l)m b 4 - 10/ (3,l,2)m c 4 - 10 h(3, 1, 2)m c 4 + 10 A(3, 2, l)m c 4 

+ 10 J (3, 2, l)m c 4 + 10 4°' 1] (3, 2, l)m b 2 + 60 3, l)m b 2 + 20 /1 0,11 (3, 2, l)m b 2 

-10 A(2, 2, 2)m fe 4 + 20 J (2, 3, l)m 6 m c 3 - 20 7 2 °' 1] (3, 2, 2)m b m c 3 + 20 f 1 (2, 2, 2)m b m c 3 

+ 10 J 2 (3, 2, l)m 6 m c 3 + 20 J 2 (2, 3, l)m 6 m c 3 + 10 A(3, 2, l)m 6 m c 3 + 20 J 2 (2, 2, 2)m b m c 3 

+20 7i(2, 3, l)m b m c 3 + 30 J (4, 1, l)m fe m c 3 + 20 /{ 0,11 (3, 2, 2)m 6 2 m c 2 + 10 J 2 (3, 2, l)m b 2 m c 2 

-20 A(3, 2, l)m 6 2 m c 2 - 30 A (4, 1, l)m b 2 m c 2 + 20 J„(3, 2, l)m fe m c 3 + 30 A (4, 1, l)m 6 m c 3 

+30 J 2 (4, 1, l)m b m c 3 + 20 J (2, 2, 2)m b m c 3 - 20 J 1 (2, 2, 2)m b 2 m c 2 - 40 J 2 (l, 3, l)m 6 m c 

-30 J 2 (2, 2, l)m b m c - 70 A(2, 2, l)m b m c - 10 7 2 0,11 (3, 2, l)m 6 m c - 20 4°' 1] (2, 2, 2)m b m c 

-20 7 2 0,11 (2, 2, 2)m fe m c + 20 7 (1, 2, 2)m fe m c + 10 7<f ' 2] (3, 2, 2)m b m c + 10 /{°' 2] (3, 2, 2)m fe m c 

-20 4 0,11 (3, 2, l)m b m c + 20 7 2 (1, 2, 2)m b m c - 20 /{°' 1] (2, 3, l)m b m c - 20 7 (2, 2, l)m b m c 

-20 7 2 0,11 (3, 1, 2)m b m c + 20 7 2 (2, 1, 2)m b m c - 40 7 (1, 3, l)m b m c - 20 (3, 1, 2)m b m c 

-20 7 2 0,1] (2, 3, l)m b m c - 10 7{ 0,1] (3, 2, l)m fe m c + 20 7 (2, 1, 2)m b m c + 10 7 2 0,2] (3, 2, 2)m fe m c 

-10 7 (1, 2, 1) - 20 7x(l, 2, 2)m 6 2 + 10 J„(2, 2, l)m fe 2 + 20 7j°' 1] (3, 1, 2)m 6 2 

-30 A (2, 1, 2)m fc 2 + 20 /J 0,11 (2, 2, 2)m 6 2 + 30 7 2 (2, 2, l)m b 2 - 10 /{°' 2] (3, 2, 2)m fc 2 

+ 10 (3, 1, 2)m c 2 - 20 7 (2, 1, 2)m c 2 - 10 4 0,11 (3, 2, l)m c 2 - 20 A(2, 1, 2)m c 2 

+60 7x(2, 2, l)m 6 2 - 30 7 2 (1, 2, 1) - 20 4°' 1] (3, 2, 2)m fe m c 3 + 20 7 1 (2, 1, 2)m b m c 

-20 4 0,1] (2, 3, l)m b m c + 20 7i(l, 2, 2)m b m c - 10 A(3, 1, l)m b m c - 20 /{°' 1] (3, 1, 2)m b m c 

-20 7{ 0,11 (2, 2, 2)m 6 m c - 40 7\(1, 3, l)m b m c - 20 /{ 0,11 (3, 2, 2)m fe m c 3 - 10 /{°' 1] (3, 2, l)m c 2 

+ 10 j} 0,1] (3, 1, 2)m c 2 - 20 7 2 (2, 2, l)m c 2 + 10 7 (3, 1, l)m c 2 + 10 7 2 0,1] (3, 2, l)m c 2 

-10 7 2 0,11 (3, 1, 2)m c 2 + 10/2(3, 1, l)m c 2 

C T - PT = -15 7 (2, 2, 2)m c 4 m fe 2 - 15 7 (4, 1, l)m c 4 m 6 2 + 5 7 (3, 1, 2)m c 5 m b - 5 7 (3, 2, l)m c 5 m b 
-5 7 (3, 2, 2)m c W + 5 7 (3, 2, 2)m c W + 5 7 (3, 2, 2)m c W - 5 7 (3, 2, 2)m c W 
+30/<f' 1] (2,2,2)m c 2 m b 2 - 15 7 (2, 1, 2)m c 2 m b 2 - 15 4° ,2] (3, 2, 2)m c 2 m fe 2 
-20 7 (2, 2, l)m 3 m b + 10 /{ 0,11 (3, 2, l)m c 3 m 6 + 10 J (2, 1, 2)m 3 m b - 30 / (1, 4, l)m c m b 5 
-5 (3, 2, 2)m c m fe 5 + 10 / (3, 2, l)m c m fe 5 + 10 J (2, 3, l)m c m 6 5 - 10 J„(3, 2, l)m c 2 m fe 4 
+30 / (1, 4, l)m c 2 m b 4 + 15 / (4, 1, l)m c 3 m 6 3 - 10 /i°' 1] (3, 2, 2)m c 3 m 6 3 
-10 / (2, 3, l)m c 3 m 6 3 + 10 / (2, 2, 2)m c 3 m b 3 + 15 4°' 1] (3, 2, 2)m c 4 m fe 2 

+ 10 / (1, 2, 2)m c m 5 3 + 15 / (2, 1, 2)m c m b 3 - 20 4°* 1] (3, 1, 2)m c m b 3 + 10 (2, 3, l)m c m b 3 
-5 / (2, 2, l)m c m 6 3 - 10 /{° ' 1] (2, 2, 2)m c m b 3 - 20 7 (1, 2, 2)m c 2 m fe 2 
+ 15 (3, 2, l)m 2 m b 2 - 5 7 (3, 1, l)m c 2 m fe 2 + 30 7 (1, 3, l)m c 2 m b 2 

+ 15/i°' 1] (3, l,2)m c 2 m 6 2 - 30 7 (2, 1, l)m c m fe + 10 4 0,1] (3, 1, l)m c m 6 - 5 4° ,2] (3, 2, l)m c m b 
-10 /i 0,11 (2, 1, 2)m c m 6 - 50 7 (1, 2, l)m c m b + 10 7 (1, 1, 2)m c m b + 5 4°' 2] (3, 1, 2)m c m fe 
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, (in ffO.ll 

+20 7<f ' 1 


(2,2,1 


—20 ' J 


(1,2,1 


-, - ?f0 21 

-15 I l ' 1 


(3,1,2 


+20 7^ ' J 


(1,2,2 


on f[0,l] 

-30 r 


(2,2,1 


-. *■ t\0 21 

-15I l ' 1 


(3,1,2 


-30 4°' 1] 


(1,4,1 


+io/r ] 


(3,2,1 



\m c m b — 70 io(l, 3, l)m c m b 3 + 5 7o°' 2 '(3, 2, 2)m c m b 3 + 5 7 (3, 1, l)m c m b 3 

- 5 4°* 21 (3, 2, l)m b 2 - 15 7 (1, 1, 2)m 6 2 + 30 4°' 1] (2, 1, 2)m b 2 
m b 2 - 30 3, l)m b 2 + 10 7 (1, 2, l)m 5 2 - 15 4°' 2] (2, 1, 2) 

— 2 I* /' ' 21 (2, 2, 2)m fe 2 + 20 7 (2, 1, l)m 5 2 + 30 7 ' 1] (2, 1, 2)m c 2 
4°* 1] (3, 1, l)m 6 2 - 5 4°' 1] (2, 2, l)m fe 2 + 20 7 (1, 2, l)m c 2 
4 0,21 (3, 2, l)m c 2 - 20 7 (1, 1, 2)m c 2 + 15 7 (2, 1, l)m c 2 



m 2 — 15 7, 
m 2 + 10 7, 



2 + 15J l [0 ' 21 



m b 4 + 15 7 (1, 1, 1) + 10 7 (2, 2, l)m fe 4 + 10 4° ' 1] (2, 2, 2)m b 4 

5 7 ° ' 2] (3, 2, 2)m b 4 + 15 7 ° ' 1] (3, 1, 2)m c 4 - 5 7 (1, 2, 2)m b 4 



-lo/„(2. 

+15 7 (2, 2, l)m c 4 + 20 /{ 0>1I (1, 1, 2) + 15 if' 21 (2, 2, 1) 
Cl~ PT = +10/ (3, l,l)m c m b + 10/l ' 1] (2,2,2)m c m 5 - 10 4°' 1] (2,2,2) 



)m r 



m c m h 



-35 7 2 (2, 2, l)m c m 6 - 10 4°' 1] (2, 2, 2)m c m 6 - 10 7 (1, 1, 2) + 20 h(l, 2, 2)m c 2 
+ 15 4°' 1] (4, 1, l)m r 2 + 25 li°' 1] (3, 1, 2)m 2 - 15 7 n (2, 1, 2)m 2 - 5 I, (3, 1, 1 W 



25 4°' 1] (3, 1, 2)m c 2 - 15 7 (2, 1, 2)m c 2 - 5 7^(3, 1, l)m c 2 
- 15 7{ 0,1] (3,1, 2)m c 2 + 30 if' 11 (2, 2, 2)m c 2 + 15 h (2, 1, 2)m c 2 + 30 4° ,1] (2, 2, 2)m c 2 
+5 7 2 °' 1] (3, 2, l)m c 2 - 15 7} ' 11 (3, 2, l)m c 2 + 15 If' 21 (3, 2, 2)m c 2 - 15 if ' 2] (3, 2, 2)m c 2 

2 



-is ir 

?[o,i]. 



>m/W- 5 7 2 (3,2,2) 



m 3 m b 3 



- 15 7{ 0,1] (3,2, l)m c - + lbir-'{ 
+ 10 7 2 (2, 2, l)m c 2 + 10 7 2 (3, 1, l)m c 2 - 5 4(3, 2, 2> 

+5 7 2 (3, 2, 2)m c 2 m fe 4 - 15 A (4, 1, l)m 3 m b + 10 7 2 (2, 2, 2)m c 3 m fe - 15 4°' 2l ( 2 , 2, 2) 
-15 7\(2, 1, 2)m c m b - 5 7 2 °' 1] (3, 2, l)m c m 6 + 5 4°' 2] (3, 2, 2)m c m b - 10 /i°' 1] (3, 2, l)m c m b 



m c m b 



-15 4 0,11 (3, 1, 2)m c m b + 15 If ' 1] (3, 1, 2)m c m b + 10 7 2 (1, 2, 2)m c m b - 10 /{°' 1] (2, 3, 1) 
+5 A(3, 1, l)m c m 5 + 30 7 2 °' 1] (2, 1, 2) - 20 7{°' 1] (2, 1, 2) + 10 7{°' 2] (3, 1, 2) 
-5 7 2 °' 1] (3, 1, 1) - 10 2, 2)m c m fe + 10 7 (1, 2, 2)m c m fe - 15 7 2 0,1] (3, 1, 2)m c m b 
-15 7 2 °' 2] (3, 1, 2) + 30 7 2 (1, 3, l)m b 2 + 20 7 2 (2, 2, l)m fe 2 + 5 7 2 (3, 1, l)m b 2 
-5 4° ' 2] (3, 2, 2)m 5 2 - 10 7 2 °' 2] (3, 2, 2)m fe 2 + 10 4 ' 1] (2, 2, 2)m 6 2 + 5 7 (2, 2, l)m b 2 
-30 4, l)m fe 2 - 10 h(2, 2, l)m 6 2 - 5 7 (1, 2 , 2)m 6 2 + 5 7{°' 2l (3, 2, 2)m fe 2 

+ 10 7 2 °' 1] (3, 2, l)m b 2 - 10 7l ' 1] (2, 2, 2)m fc 2 - 15 7 2 (1, 2 , 2)m b 2 + 30 7{°' 1] (1, 4, l)m 6 2 
-5 7?' 1] (2, 2, 1) - 5 7 2 °' 1] (2, 2, 1) + 20 7«f ' 1] (1, 2, 2) - 15 7 2 (2, 1, 2)m 6 2 - 30 7 2 °' 1] (1, 4, l)m 6 2 
+ 15 7 2 °' 1] (3, 1, 2)m b 2 + 20 7 2 °' 1] (2, 2, 2)m fe 2 - 10 7{°' 1] (3, 2, l)m b 2 + 15 If ' 1] (3, 2, l)m fe 2 
+5 7!(l,2,2)m fe 2 + 15 7}°' 21 (2, 2, 2) - 15 7 2 0,2] (2, 2, 2) + 10 4°' 1] (3, 1, 1) 
+ 15 (2, 2, 1) - 5 7{ 0,3] (3, 2, 2) + 5 7 2 M (3, 2, 2) - 5 7 2 (3, 2, 2)m c 6 - 5 7 (3, 2, 2) 
+5 A(3, 2, l)m c 4 - 10 7 (3, 2, l)m c 4 - 15 7 2 (4, 1, l)m c 4 - 15 7 2 (2, 2, 2)m L ' 



-15 7!(2,2,2)m 



15 A (4, 1, l)m c 4 + 15 7 2 °' 1] (3, 2, 2)m c 4 - 15 7{ 0,1] (3, 2, 2)m c 



-15 7 (2, 2, 2)m c 4 + 5 7 2 (3, 2, 2)m c 5 m b + 5 7 (3, 2, 2)m c 5 m fe - 5 A(3, 2, 2)m c 5 m b 
+5 7 (3, 2, 2)m c 4 m fe 2 - 5 7^3, 2, 2)m c 4 m fe 2 + 5 7^3, 2, 2)m c 3 m fe 3 + 5 7^3, 2, 2)m c 6 
-5 7i(3, 1, 2)m 3 m b + 10 If' 11 (3, 2, 2)m c 3 m 6 + 10 A(2, 3, l)m 3 m b - 10 4°' 1] (3, 2, 2)m c 3 m fe 
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+5 7 (3, 1, 2)m c 3 m b - 5 A(3, 2, l)m c 3 m fe - 10 h(2, 2, 2)m c 3 m b + 5 7 2 (3, 1, 2)m c 3 m b 
+ 15 7 2 (4, 1, l)m c 3 m 6 + 10 7 (3, 2, l)m c 3 m 6 - 10 7 2 0,11 (3, 2, 2)m c 3 m b + 5 j|°' 1] (3, 2, 2)m 6 4 
-15 7 (4, 1, l)m c 4 + 5 7\(3, 1, 2)m c 4 - 10 J 2 (3, 1, 2)m c 4 - 5 J (3, 1, 2)m c 4 
J (2, 2, 2)m 3 m b + 5 J 2 (3, 2, l)m c 3 m b - 10 J 2 (2, 3, l)m c 3 m b - 10 J 2 (2, 2, 2)m c 2 m b 2 
-20 J 2 (3, 2, l)m c 2 m 6 2 + 30 J (l, 4, l)m c 2 m b 2 + 10 h{3, 2, l)m 2 m b 2 - 5 J (3, 2, l)m c 2 m 6 2 
+ 10 (3, 2, 2)m c 2 m b 2 + 30 7 2 (1, 4, l)m fc 4 - 15 J 2 (4, 1, l)m 2 m b 2 + 30 7 2 (1, 4, l)m c 2 m fc 2 
+5 7 2 (3, 2, l)m c m b 3 + 5 7{°' 1] (3, 2, 2)m c m b 3 + 5 7 2 (3, 1, 2)m c m b 3 - 30 7 2 (1, 4, l)m c m b 3 
-10 7 2 (2, 3, l)m c m b 3 - 5 7 2 (3, 2, l)m b 4 - 30 7 (1, 4, l)m c m 6 3 + 10 7 (2 , 3, l)m c m b 3 
+ 10 7 2 °' 1] (2, 3, l)m c m b - 40 7 (1, 3, l)m c m b + 15 7 (2, 1, 2)m c m b - 10 7 2 (2, 2, 2)m 6 4 

where 

(M 2 ) 1 (M 2 ) 3 I n (a,b,c) . 



a,b,c) = (M 2 



Ml 



a" 



d ( M?Y d. ( M?V 
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Figure 7: The decay width for B c — > D sl {x + {i~ with respect to 9 S and q 2 . The left 
figure shows decay width of B c — > _D s i(2460)/i + /i~ and right figure belongs to B c — > 
D s i(2536)// + // _ . In these figures, x = 9 S , y = q 2 and z = r(S c — > D s ifi + fi~) x 10~ 20 . 




Figure 8: The decay width for B c — > D s it + t with respect to # s and g 2 . The left 
figure shows decay width of £? c — > £> s i(2460)t + t~ and right figure belongs to S c — > 
D s1 (2536)t + t-. In these figures, x = 6 S , y = q 2 and z = r(£ c -> D s1 t + t~) x 10 -21 . 
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Figure 11: The same as Figs. 9, but for B c — > D s \vv. 
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